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In this paper we study the problem of wave propagation in a 3-D optical fiber. The goal 
is to obtain a solution for the time-harmonic field caused by a source in a cylindrically 
symmetric waveguide. The geometry of the problem, corresponding to an open waveg- 
uide, makes the problem challenging. To solve it, we construct a transform theory which 
is a nontrivial generalization of a method for solving a 2-D version of this problem given 
by Magnanini and Santosa0 

The extension to 3-D is made complicated by the fact that the resulting eigenvalue 
problem defining the transform kernel is singular both at the origin and at infinity. The 
singularities require the investigation of the behavior of the solutions of the eigenvalue 
problem. Moreover, the derivation of the transform formulas needed to solve the wave 
propagation problem involves nontrivial calculations. 

The paper provides a complete description on how to construct the solution to the 
wave propagation problem in a 3-D optical waveguide with cylindrical symmetry. A 
follow-up article will study the particular cases of a step-index fiber and of a coaxial 
waveguide. In those cases we will obtain concrete formulas for the field and numerical 
examples. 

Keywords: Wave propagation. Optical waveguides, Helmholtz equation. Green's func- 
tion. Spectral representation. 

1. Introduction 

In this paper, we study the wave propagation in a cyhndrical optical fiber." As 
model equation, we use the Helmholtz equation 



also called the time-harmonic wave equation. The number k is called the wavenum- 
ber, and the function / represents a source of energy. We require that the index of 
refraction n{x, y, z) have the form: 



where R is the radius of the waveguide. 

The main result of this paper is the construction of a representation formula for 
a solution u of satisfying suitable radiation conditions. Our results generalize 
a similar formula obtained by Magnanini and Santos in the two dimensional case. 

In El it is shown that the energy of the electromagnetic field is divided into two 
parts: a part propagates inside the waveguide as a finite number of distinct guided 
modes, while the other part either decays exponentially along the fiber or is radiated 
outside. Our case reveals a new feature: for special choices of the parameters, new 

"^We will use the terms optical waveguide and optical fiber interchangeably. 



Au + k^n{x, y, z)'^u 



f{x,y,z), {x,y,z)eR^, 



(1.1) 
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kinds of guided modes appear which, rather than decaying exponentiaUy outside 
the fiber, vanish as a power of the distance from the fiber's axis. 

As inj^ we use the technique of separation of variables. But there are a number 
of important differences. We separate variables in the three cylindrical coordinates 
r, z9, and z. The component in is simple to solve for, and will only introduce 
a sum over m G Z. The separation in z can be solved as in.l^ More complicated 
is the study of the r-coordinate: we will obtain a differential equation in r which 
is Bessel's equation, save for a non-constant coefficient. This equation will have a 
singularity at r = 0. Also, instead of working as in 2-D with sine and cosine functions 
whose zeros are uniformly distributed, we will work with Bessel's functions whose 
distribution of zeros is more complicated. These differences make the method used 
in 1^ inapplicable. 

Instead, we use the theory of singular self-adjoint eigenvalue problems for second 
order differential equations as presented in El and. ^ Applying this theory to our 
problem and doing the explicit calculations took some effort, chiefly for the reason 
that our self-adjoint eigenvalue problem has (just like in the 2-D case) a coefficient 
which is, with some restrictions, a general function, so we cannot obtain solutions 
for our equation in terms of concrete functions. 

The paper is organized as follows. In section 2 we derive the second order 
self-adjoint eigenvalue problem associated with the Helmholtz equation H1.I|I . In 
section 3 we will prove a set of technical lemmas aimed at studying the behavior of 
the solutions of this eigenvalue problem. In section 4 we will classify the solutions of 
the eigenvalue problem which are "well-behaved" (in a sense to be specified there) 
as r — > and r oo. The motivation is that the electromagnetic field in the fiber 
will have a representation in terms of the "well-behaved" solutions of this eigenvalue 
problem. In section 5 we summarize the theory of self-adjoint eigenvalue problems 
as exposed inl^and.l^ 

The functions defined in section 5 are calculated in section 6. In section 7 the 
transform defined in section 5 is computed. The obtained transform is used in 
section 8 to find Green's function for the Helmholtz equation 1)1. l|l , and in turn, to 
find the desired electromagnetic field in the fiber given the source. 

2. The Eigenvalue Problem 

A typical optical fiber is a cylindrical dielectric waveguide, made of silica glass or 
plastic. Its central region is called core, surrounded by cladding, which has a slightly 
lower index of refraction. The cladding is surrounded by a protective jacket. Most 
of the electromagnetic radiation propagates along the core. The electromagnetic 
field intensity in the cladding decays exponentially along the radial direction. This 
is why, the radius of the cladding, which is typically several times larger than the 
radius of the core, can be considered infinite. 

A typical optical fiber is weakly guided, which means that the difference in the 
indices of refraction of its the core and cladding is very small. In these conditions the 
electromagnetic field in the fiber is essentially transverse with each of the transverse 
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Figure 1: An optical fiber. 

components approximately satisfying the Helmholtz equation Hl.l|) . This is the 
so-called weakly guided approximation. For a derivation of Hl.ll) from Maxwell's 
equations for a weakly guided fiber seeP chapters 30 and 32. 

Because of the cylindrical geometry it will be convenient to use the cylindrical 
coordinate system (r, z), with z being the axial direction. Then, the index of 
refraction will depend on the r variable only, n = n{r). In the new coordinates 11.11) 
becomes 

d'^u 1 d { du\ 1 d'^u ,9 , „ N /„ N 

Jz^^rd-rV^rj^r^W^^ ^^^^ = ^'''^ 

This is a linear partial differential equation. The solution of this equation will be 
determined as soon as we find its Green's function. In order to obtain the latter, 
consider the homogeneous version of H2.1fl . 

d'^u I d / du\ 1 d'^u ^2 / \2 Q 
dz^ r dr\ dr ) d'd'^ 

Look for a solution in separated variables, u(r, t9, z) = Z{z) v{r). It is quickly 
found that we must have 

it(r,t?,z) = e*^'=^e""''z;(r-), (2.2) 
with P (z C, m ^ Z, and v{r) satisfying the differential equation 

v" + -v' + ikMrf ~l3^-^\v = 
r I r^ } 

(the derivative here, and in the rest of this paper will always be in respect to the 
variable r). 

Let i? > be the radius of the fiber core. Then n{r) = nd for r > R. Denote 
d'=nl-n% l = e{nl-P^), q{r) = e[nl - n{rf]. (2.3) 
Then this equation becomes 

If 2 ^ 

v" + --v' + U-q{r)- — \v^O. (2.4) 
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We will view (|2.4(l as an eigenvalue problem in ^ e C. The variable r is in (0, oo), 
the number m is an integer, the function q is bounded, measurable, real- valued and 
non-negative, with q(r) — (P > for r > R > 0. It will be convenient to make a 
variable change. Denote w(r) — y/rv{r). Get the equation 

w" + (^l-qir)-^^^--^^w^O, re(0,oo). (2.5) 

This will be our self-adjoint eigenvalue problem. 




Figure 2: The function q(r). 



3. A Study of the Solutions of the Eigenvalue Problem 

Before going further we will need some information about the behavior of the 
solutions of the differential equation (|2.5II as functions of r and I and m. This will 
be the subject of the next four lemmas. 

Lemma 3.1 There exists a solution jm(r, ^) (r > 0, ^ € C, m e Z) of (|2.5|l sucii 
that 

hm - 1. lim „ , ^'"/'"'^} , =1. (3.1) 

r^0rl™l+V2 ' r^o (|m| + 1/2) rl™l-i/2 ^ ' 

The functions jmir, I) andj^(r,Z) are analytic in I asr is fixed. There exists another 
solution yjn{r,l) of (|2.5|l such that 

lim ^^pM^ = 1, lim = 1, if \m\ > 1, (3.2a) 

and 

lim^I^iM^l, hm /^/,^'^^, =1, ifm^O. (3.2b) 

r^o -^rmr i-^o lnr/(2Y'') 

Proof. We will assume to > 0, and then set j^m = jm and y-m = y-m- Make the 
variable change w — r^+^/^cr in ||2.5|) . Obtain 

2m + 1 , , , ^ 

cr" —(t' + {l- q{r)} cr = 0. 

r 

Denote k = 2m -t- 1, fc > 1. Multiply this equation by r*^. Get 

{r'^a'Y = r''iq{r)-l)a. (3.3) 
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To prove this lemma we need to find two solutions a{r, I) and r(r, I) of H3.8|l such 
that 

limcr(r,0 = l, limcr'(r,/) = 0, (3.4) 



mil r^™-''- '"1 — 1 ^2m+l^/ 

and 



lim r2"T(r, = 1, lim r^"'+W'{r, I) = -2m, if m > 1, (3.5a) 



lim = 1, lim rr' (r, I) = 1, if m = 0. (3.5b) 

r— Inr r^O 

Also, (7(r, ^) must be analytic in I for r fixed. The idea to finding a and r is to 
rewrite (13.3(1 as an integral equation. 

Let = sup,^g[Q Q^-j q{r). If lu : [0, oo) ^ C is a function, bounded and integrable 
on every compact subset of [0, oo), and p > is an integer, then one has 

r r 



sPuj{s)ds < sup \uj(t)\ sPds^ sup \uj{t)\ -. (3.6) 

te[o,r] J te[0.r] J3 + 1 

u 

Consider the operator 

r t 

Tuj{r)= Jt-'' Js''{q{s)-l)cj{s)dsdt, (3.7) 



defined for complex-valued functions lu which are bounded and integrable on every 
compact subset of [0,oo). By applying H3.6|l to the inner-most integral in H3.7|l we 
deduce 



sup \TLu{t)\ < (goo + 1^1) sup \Lu{t)\- 



r2 



te[0,r] t6[0,r] 2 (fc -f 1) 

In particular, for any r > 0, T is a bounded linear operator from the space C([0, r]) 
of continuous, complex-valued functions defined on [0, r] onto itself. By using the 
same reasoning as above, one can show by induction that for any integer n > 

|T"^(r)| <(goo + |^|)" sup \Luit)\- 



t^lnr] " {2 ■ 4 - ■ - 271} {{k + 1) ■ {k + 3) ■ ■ ■ {k + 2n - I)}' 

(3.8) 

Using (|3.fi|l it is easy to check that if function a satisfies (|3.4(l . then 1)3. 3|l is 
equivalent to 

r 

<j'{r, A) = 1 / /-(g(.s) - A) a{s, A) ds, (3.9) 



which in turn is the same as 

a = l + Ta. (3.10) 

We will find a by applying the method of successive approximations. Let ao = 1, 
and cr„+i = (To + Ta„, n > 0. Then, 

(T„ = (To + Tao + ■ • • + T"cto, n>0. 
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By using 1)3. 8|l we obtain that the series 

a = (Ta+Tao + --- + T"(Ta + ^0 = 1, (3.11) 

is uniformly convergent for r and I in compact sets. Since f7„ is continuous in r and 
analytic in I, the same property will hold for a. Using the fact that the operator T 
is continuous on C([0, r]) for any r > 0, it follows from the above series that 1)3.10(1 
holds, and therefore (|3.3|l holds. 

Now prove the existence of r satisfying (|3.3|l with the boundary conditions (|3.5a(l 
or ((3.5b|l . That will be done in several steps to be described below. For the rest 
of the lemma, a = a{r, I) will be the solution to 1)3.3(1 found above. We are not 
looking for any properties of r in respect with I. Therefore, we will fix Z G C, and 
will consider r to be a function of r only. At step 1 we will prove that if t is a 
solution to 1)3.3(1 which is linearly independent of a, then t is unbounded as r — > 0. 
At step 2, we will show that any solution r to ()3.3(l has the property 

\T{r)\ < A + Br-'',0 <r <ro, (3.12) 

for some ^>0,i3>0,ro>0. At step 3 we will prove that any solution r of 1(3.3(1 
linearly independent of a satisfies 

r'(r) = 0{r-'') as r ^ 0. (3.13) 

That will imply l)3.5a(l or ()3.5b(l depending on whether to > 1 or 771 = 0, that is, 
A: > 3 or fc = 1. 

Step 1. Show that any solution t to 1(3.3(1 which is linearly independent of a is 
unbounded as r ^ 0. Assume that t is such a solution and that it is bounded as 
r — > 0. From 1(3.3(1 we get 

r 

r^T'{r) = ci + j s^{q{s) - l)T{r) ds 


for some constant ci. We cannot have ci ^ 0, because then T'(r) = 0(r^'^) with 
k > 1, so r cannot be bounded as r ^ 0. Then, if Ci = we can divide by r'^ and 
integrate again, to obtain 

r = C2 + Tt, 

with T the operator defined by 13.7(1 and C2 another constant. Since the solution a 
satisfies a = 1 + Ta, we deduce that — t — cio will satisfy 

Then ip = T"(^ for any ri > 0. By using ((3.8(1 we get 

te[o,r] {2-4---2n|{(fc+l)- (fc + 3) • • • (fc + 2n - 1)} 
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for any r > and n > 0, which imphes ip = 0, that is, r — C20. This is a contradiction 
with the assumption that r is hnearly independent of a. 

Step 2. Show that for any r solution of inequahty l|8.12(l holds. Let us 

again write 1)3. 3|l as an integral equation. This time we cannot integrate from r = 0, 
since we expect an unbounded solution as r ^ 0. Let 7'o > be a fixed number. It 
is easy to show that r will satisfy 1)3.3(1 if and only if 



r = To + St, 



(3.14) 



where 



T^ir) ^co + do It ''dt 



with Co = ''"(^0)1 do ~ rQr'(ro), and the operator S is defined on functions 9 inte- 
grable on compact subsets of (0, 00) and is given by the formula 



Se{r)= t-'' s''{q{s) -l)e{s)dsdt. 



Use once again the method of successive approximations to express t as the 
sum of an infinite series. Let co G M, do G and define To(r) = co + rfo J^t-'^dt, 
Tn+i = To + Stji, n > 0. We have 

T„ = To + S'to + • • • + 5"to, 71 > 0. 

Notice that if uj is an integrable function on (0, tq), then for < r < ro 



s ''Lu{s)ds 



< I s-''\uj{s)\ds Kr-" l\uj{s)\ds. 



-k 



(3.15) 



By applying <|3.15|l we get the estimate |To(r)| < |co| + |do|''o?' for < r < rg. 
Then, 



fc 



s'=(g(s) - /)to(s) ds 



< /('Zoo + |^|)(r^1co| + Molr-o) ds = C{q^ + \l\){ro - t), 



where C = Tq |co| + |do|''o- Apply l|3.15|l to estimate S'to. 



\SMr)\ < r-^- / C{qoo + \l\){ro ^ t) dt = C{q^ + \l\)r 



2! 



By repeating the above reasoning, one obtains by induction 

|^"ro(r)| < C(qoo + \l\Tr-^ ^ ° , , < r < ro, n > 1. 

(2n)! 
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As it was done for cr, one can prove that the series 



E 

Tl = 



is uniformly convergent for r in compact subsets of (0, vq] and that t satisfies H3.14(l . 
One can then estimate r. 



oo 



r(r)| < y I^^ToMI < |co| + Cr-*=eV9-+l'K'-o-'^), < r < ro, 



n=0 



which imphes l|3.12(l . 

Step 3. Show that any solution r to (|3.3(l which is linearly independent of a 
satisfies (jH.lHII . From H3.12|l we deduce that r'^{q{r) — l)T{r) is bounded as r — > 0. 
Integrate H3.3|) from to r. Get 

r 

rV(r) - C3 + / s'^iqis) ~ l)Tis) ds, (3.16) 



with C3 a constant. 

Suppose that C3 = 0. Then, from (pDUl and (pTT^ it follows that |r'=r'(r)| < Mr 
for < r < ro and some M > 0. But then, |T'(r)| < M/r''"^ for < r < rg. Note 
that fc is an odd number, since k = 2m + 1 with m > an integer. In particular 
k 2, OY k — 1 ^ 1. Since r(r) = r(ro) + J^^t'{s) ds, we obtain 

Kr)|<Ai + -^, 

for some Ai > 0, Bi > 0. This is the same as (|3.12ll but with k — 2 instead of k. By 
repeating several times the same reasoning starting with (|3.16() and the assumption 
C3 = we will be able to reduce the exponent of r in this inequality by 2 every 
time, until we get that t must be bounded as r 0. As shown at step 1, then r 
is linearly dependent of a, which is a contradiction. Thus C3 7^ 0. Then, (|3.12|) and 
(pnH|l give us r'(r) = 0{r-'') as r ^ 0, or ^J^. 

From H3.13|l and the equality T(r) = T(ro) + J^^t'{s) ds it follows that as r ^ 0, 
T(r) — O(lnr) if = 1, and r(r) — 0{r^'^^^) if fc > 1. One can then get either 
(|3.5a|) or Ij3.5b|) by conveniently multiplying r by a non-zero constant. This proves 
the lemma. □ 

The following lemma will study the properties of jmif, ^'^^ ^ real. 

Lemma 3.2 Let / = A e M. Then jm{r,X) and j^(r, A) are reai. Jf A < 0, then 
j„(r, A) > rl™l+V2 and j'^{r, A) > 0. 

Proof. We will consider only the case m > 0, since jmir, A) is an even function of 
m. If A is real, then q{r) — A is real. With the notation from the proof of lemma 
13.11 jmir, A) = r"^^^^^a{r, A), with a{r, A) given by the series (|3.11() . By using this 
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series, and definition H3.7|l of the operator T we deduce that a{r, A) is real, and 
consequently, jmii", A) is real. 

If A < 0, then q{r) — A > 0. In this case, the operator T will map non-negative 
functions into non-negative functions. We conclude from the construction (13.11(1 of 
cr that cr(r, A) > 1, and so, j,n{r, A) > r"+i/2. 

To show A) > it suffices to prove a'(r,X) > 0, which is an immediate 

consequence of ((3.9(1 . □ 

The next lemma will study the properties of j„t(''iO foi' large absolute value of 
m and I — X a real number. 

Lemma 3.3 Let r > be a fixed number and A C [0, oo) be a bounded set. If \m\ 
is sufEciently large, then for all X E A 

j™(r,A) > 0, j'^ir,X) > 0. 

Proof. We can assume to > 0. Denote k — 2m + 1 > 1. We have that jmir, A) — 
r™+^/^(T(r, A), with cr(r, A) defined by the series 1(3.11(1 . Here we will prefer to use 
the notation crm(r. A) to emphasize its dependence on to. By using this series and 
estimate ((3.8(1 it follows that 

|1 - fT™(r, A)| < AeA, 

with C = C(A, qoc) a constant, where goo = sup^gjQ o^-) q{r). 
From 113.9(1 we deduce 

WLirA)\ < AG A, 

with D — D(A, Qoo) another constant. 

These observations imply that for to large enough and A G A, am{r,X) is close 
to 1, while (y'jn(r, A) is close to zero. Then, since 

j-^(r,A)=r™+i/V„(r, A), 

and 

j'^ir, A) - r™-i/2{(^ + 1/2) a,„(r, A) +ra'„,ir, A)}, 

it follows quickly that for to large enough both of these quantities are strictly 
positive. □ 

Lemma 3.4 For I E {—tt/2 < arg(Z — d^) < 3tt/2}, there exist two solutions 
Wmi'T, I) and Xmir, I), r > m E Z, of ((2.5(1 which are analytic in I. If I E C+ = 
{I : Jml > 0}, then as r — > oo, Wm{r,l) together with its derivative will decay 
exponentially, while Xm{r, I) and its derivative will increase exponentially. 

Proof. Recall that q{r) = for r > R, therefore on [i?, oo) 1(2.5(1 becomes 
w" + U-d^ U = 0. 
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Set 

wUr,l)^V^H^^'>{Vl~d^r), x^rj) ^ V^H^^^l ~ r), (3.17) 

where H^\C) andi/i''(C) are the Hankel functions of m-th order. They are analytic 
functions with the domain {— tt < argC < tt}. Consider the function y/z defined on 
{— 7r/2 < argz < 3tt/2} with values in {—tt/A < argC < 37r/4}. We obtain that 
Wm{r,l) and Xjn{r,l) are analytic functions of / £ {—tt/2 < arg(Z — (P) < 37r/2}, 
which is the whole complex plane except those I for which I ~ d? has a zero real 
part and a non-positive imaginary part. 

According to formulas (9.2.3) and (9.2.4) from;IJwe have the asymptotic expan- 
sions 

H,W(C) - V27K)e'^'^"'""^'""/^^ 

and 

H^i\0 - V27R)e-^(^-""/'-"/') 

as ICI oo and | argCI < tt. Then, 

u;„(r,0 = O(e'^^0' xUr,l)^0{e-'^^'') as r ^ oo. 

If I G C+, then Jm — cP > and thus, i^/l — (P has a strictly negative real 
part. Therefore, as r — > oo, w,„(r, Z) will decay exponentially, while Xmir^l) will 
increase exponentially. By formally differentiating the above equalities (for a rigor- 
ous justification one needs to use equalities (9.2.13) and (9.2.14) from^) we deduce 
that w'^^^{r,l) will decay exponentially, while x'^^{r,l) will increase exponentially as 
r — > oo. □ 



4. Classification of the Solutions 

The main purpose of this paper is to prove that under certain conditions, the 
solution to the Helmholtz equation H2.1|l is a superposition of functions of the form 
H2.2|l . For each of the functions in the superposition, v{r) will satisfy (|2.4(1 with the 
notations of (|2.3I) . and I will be a real variable, which we will denote by A (thus 
we will reserve the notation I for the complex variable, and the notation A for its 
restriction to the real axis). In addition the following properties will hold: 

v{r) is bounded as r ^ 0, (4.1a) 

and 

_ / is in i2(i^^oo), if A < 

r y/rv[r) < (4-lb) 
I is bounded as r ^ oo, if A > d . 

In this section we will study and classify the functions v{r) with the properties 
(|4.1a|) and Ij4.1b|) . 

A solution v{r) of H2.4|) has the form v{r) — ■w{r)/y/r, with w{r) satisfying (|2.5|) . 
Lemma [3.11 shows how the solutions of (I2.5|l look like. It is clear that in order that 
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v{r) satisfy l|4.1a|l . we need w{r) = jjn{r,X), or 

y(^r) = tl^^. (4.2) 
Condition Ij4.1b(l is then satisfied if and only if 

I IS bounded as r — > oo, it A > a . 

Next we will investigate for which A condition (|4.3|l holds. 

We will need to consider four cases: A < 0, < A < d^, A = and A > d^. In 
each of these intervals jm(?', A) and consequently v{r), will have a different behavior. 

Case 1. If A < 0, then according to lemma jm{r,\) > 7-l™l+i/^, therefore 
jmii", A) will be not square integrable on (i?, oo). 

Case 2. Assume < A < c?^. For r > R the function q{r) defined in 1)2. 3|l is 
constant and equal to (P. Then 1)2. 5|l becomes 

A-d2_!^^!_iZl\u> = 0, re[i?,oo). (4.4) 



The solutions to this equation are 

kmir, A) = V^i^™(Vrf2 - Ar), \<d\r>R, (4.5) 

and 

y/PlmiVd^ - \r), \<d^,r>R, 



where Km and /,„ are the modified Bessel functions. Formulas (9.7.1) and (9.7.2) 

isions 

^/sKm{s) ^ ■\/7r/2 as s ^ oo, 



from ^ give us the expansions 



and 

\/sIm[s) ~ y/iTple^ as s — ^ oo. 

Thus, we have one solution decaying exponentially while another increasing 
exponentially. In order that jm{r, A) be in L^{R, oo) we need 

jmir, A) = Ckm{r, A) for r > R, 

for some constant C. S incc both j^^j(r, A) and Ckmij"^^ satisfy the same second 
order differential equation, namely (|4.4|) . to ask for the equality of these functions 
is the same as to ask for them to satisfy the same boundary conditions at r = i?, 
which translates into 

j™(i?,A)-Cfc™(i?,A) = 0, 

and 

j,:„(i?,A)-Cfc;,(i?,A) = 0. 
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Thus, C must satisfy simultaneously two different conditions. This is possible if 
and only if 



Then we will have 

^™("'^) = r4i4T^"("'^)' (4.7) 

krn(R, A) 

It will be shown later that for each m, the set of A such that (|4.6|) holds is finite. 
Notice that in this case jm(r, A) will decay exponentially as r ^ oo. Its derivative 
has the same property, this follows from the asymptotic expansion 

[y/sKm{s)]' y/TT/2e^'' as s ^ oo (4.8) 

(according to the formulas (9.7.2) and (9.7.4) from.-*^) In conclusion, the only A e 
(0, d^) for which (|4.3|) holds, are those satisfying H4.6|l . 

Case 3. Let now X = (P. Two linear independent solutions of (|4.4|) are in this 
case 

and 

{y^lnr for m — 0, 
^l/2+|m| fQj. ^jLQ_ 

The second solution is not bounded for any m G Z. By matching the boundary 
conditions at r = i? as in the previous case, it is easy to show that jmif^ A) will be 
proportional to the first of these two solutions if and only if 

Jm[U,l) H 

and then, 

The function jm{r, A) will be in L^{R, oo) if and only if |m| > 2. 

A function of the form H2.2|) . with v{r) given by 1)4. 2|l . for which < A < 
and either 1)4.6(1 . or 1)4.9(1 (with \m\ > 2) holds, is called a guided mode. Note that 
a guided mode decays in r either exponentially, or as r^'™' (|m| > 2). 

Case 4. Let \> d^. Two solutions of ()4.4)l are then 



a™(r,A) = V^Jm{V\~d^ r) , 5™ (r. A) = ^^ 1™ ( VA^ r) , A > d^, (4.11) 

where Jm and Ym are the Bessel functions of the first and second kind. Note the 
formulas 

^Jm{s) = V'7i72cos(s - TO7r/2 - 7r/4) + 0(s"^/2), 
\/ir„(.s) = V7i72sin(s - TO7r/2 - 7r/4) + 0(s"^/^), 
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as s — > oo (they are a particular case of formulas (9.2.1) and (9.2.2) from.-'^) We infer 
that am(r, A) and 6,„(r, A) will be bounded as r ^ oo. By formally differentiating 
the above formulas it follows that a'„(r, A) and 6^(r, A) will also be bounded as 
r — > oo. 

The functions am(r, A) and hm{r,\) are linearly independent, since Jm and 1^ 
are linearly independent. Then, for r > R, jmif, A) will be a linear combination of 
them. To find the coefficients of the linear combination, set 

j„i(r, A) = c,„(A)a,„(r, A) + (i„j(A)6,„(r, A), r > R. 

By matching the boundary conditions at r = i?, we obtain a linear system which 
enables us to solve for Cm and dm- Apply the equality 

Y:^{z)Jm{z) ~ J'rA^)Ym{z) = 2/(^z), 

((9.1.16) from ^ to find a value for the determinant of the this linear system. 
Obtain 

6'„(i?,A)a™(i?,A)-a:„(ii',A)6™(ii',A) = 2/7r, (4.12) 

and therefore, 

c™(A) = |{6'„,(i?,A)j™(i?, A) - j:„(i?,A)&„(i?,A)}, (4.13a) 

rfm(A) = -|{a:„(i?, A)i„(i?, A) - j:„(i?,A)a„,(i?,A)}. (4.13b) 

It is easy to see that jmif, A) and its derivative will be bounded as r — > oo, and 
so, condition H4.3I) will be satisfied for all \> cP . 

Let us look at the expression of (|2.2|) for \ > dP and with v{r) given by (|4.2I) . 
Notice that if < A < k'^n'^, then 12.211 will be oscillatory in z (recall that = 
fc^riQ — A). In this case we will say that (|2.2|) is a radiation mode. On the other 
hand, if A > fc^riQ, then 13 becomes imaginary. Depending on the sign of Jm/3 we will 
have exponential decay in one of the directions z ^ — oo, z — > oo, and exponential 
growth in the other one. For A > fc^riQ, (|2.2() will be called an evanescent mode. 

5. The Theory of Eigenvalue Problems 

Consider the eigenvalue problem 

w" + {l-Q{r)}w ^Q, re(0,oo), (5.1) 

where / G C. Assume that Q is integrable over any compact subset of (0, oo) (in El 
and El the theory is developed only for continuous functions Q, but it is mentioned 
in a footnote on page 224 of El that it suffices for Q to be as we assume above). 
Let < i? < oo be an arbitrary but fixed number. Let i.p{r, I) and 9{r, I) be the 
solutions of (|5.1() with the boundary conditions 

('^(i?,0=0, ^'(i?,0--l, ,J.2^ 

\e{R,i) = 1, 6i'(i?,o = 0. ^ ■ ' 
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Since (|5.1|l has an analytic dependence on the parameter I, the solutions d{r, I) and 
f{r,l) will be analytic functions of / for r fixed. 

Any solution to l|5.1|l linearly independent of ip can be represented, up to a 
constant multiple, in the form 

^P = e + Mip, (5.3) 

with M G C. Let r e M and < t < oo. Look for a solution of H5.1|l of the form 
(|5.3(l to satisfy the boundary condition 

cosTip{t, I) + sinTt/j^t, I) = 0. 
It is a direct calculation to check that we need 

M = Mil) = _^(M)cosr + ^'(M)sin.^ 

^p{t, I) cos T + if' (t, I) SlUT 

Let C+ be the open upper complex-half-plane, C+ — {I : 3ml > 0}. As shown 
in chapter 2 the following results hold: as t ^ 0, M{1) converges uniformly on 
compact subsets of C+ to a function Mo{l) analytic on C+. Moreover, the function 

MrJ)^Oir,l) + Mo{l)p{rJ), leC+, (5.5) 

is in L^{0, R), and one has 



R 



lipoirjjl dr = — — - — . (5.6) 
Jmt 

As t oo, M{1) converges uniformly on compact subsets of C+ to an analytic 
function Moo{l) on C+, and if 

^ooirJ)^9ir,l) + Moo{l)ip{r,l), leC+, (5.7) 

then tpoo{r,l) € L^(i?, oo) and 

2. 3m Moo (0 

|V'oo(r,OI dr = — . (5.8) 

Jm I 

Note that the obtained Mq, Mqo, ipo, ipoo depend on the parameter r G M. Thus, 
possibly these quantities, and therefore the representation given below, in Theorem 
15.11 will not be unique. This might be true in general, but in our concrete case, 
given by (|2.5|l . these quantities will turn out to be unique, as we will see from lemma 
16.11 So then the transform we are looking for (which is calculated in Theorems 17.11 

and l7.2ll will be unique. 

9 '7' 

In section 9.5 of and chapter 3 of'-'-' it is proved that for any A G M the following 
limits exist 

A 

C(A) = lim / -3m ds, (5.9a) 

5^0+7 Mn(s + idj~Moc(s + io) 
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A 

r,{X)= lim /-Jm „^, ^^tfl ^^ds, (5.9b) 




A 



5^0+7 A/o(s + i(5)-Moo(s + «(S) 







It is shown there that the functions ^ and C are non-decreasing, and that r] is with 
bounded variation. In addition, for any Aq < Ai real numbers, 

M^i) - ^(Ao)}' < {^(Ai) - e(Ao)}{C(Ai) - C(Ao)}, (5.10) 

as stated on page 252 of El (with a different notation). Then, equaUties (3.1.8), 

(3.1.9), (3.1.10) from^lgive us an expansion formula for a function g € L^{0, oo) in 

terms of 9{r, I), ip{r, I) and the functions ^, rj and C,. The same result is proved in El 

at Theorem 5.2. In this reference the representation result is stated more rigorously, 
7 

so we will prefer it over. To state the result we need some notations. 

Denote p = (C, ?/, C)- For any vector T = (ri,r2), where Fi, Tj : M ^ C, let 

OO 

l|r||'= y |ri(A)|2de + 2Mri(A)f2(A)}d77 + |r2(A)|2dc. (5.ii) 

— oo 

The fact that and 77 are non-decreasing, together with (|5.10l) . gives us that | |r| p > 
0. It is easy to check that || • || is a semi-norm. Denote by L^{p) the space of all 
r = (ri,r2) such that ||r|| < oo. This is then the statement of Theorem 5.2 fromPI 

Theorem 5.1 If g e L^(0,(X)), the vector T = (ri,r2), wiiere 

CSO 00 

ri(A) = J 9{r, X)g{r) dr, T^iX) = j ^(r, X)g{r) dr, 

converges in L^{p), that is, there exists T e L^{p) such that 

||r-r='*|| ^0 asc^O, d^oo, 

where for < c < d < 00 

d d 

rf (A) = j e{r, A)g(r) dr, r^'^(A) = J ^{r, \)g{r) dr. (5.12) 

c c 

The expansion 

00 

9{r) = - / {d{r, A)ri(A) ^^(A) + 0(r, A)r2(A) dr^{X) 



TT 



^(r, A)ri (A) dTj{X) + ^{r, X)T2{X) dC{X)} 
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holds, with the latter integral convergent in L^{0, oo), that is, g°"^ g in L^{0, oo) 
as cr — !■ — oo, T — > oo, where for — oo < t < a < oo 



- ^ / ^)ri(A) daX) + 0{r, A)r2(A) dry(A) 

(T 

+ ^{r, A)ri(A) dr,(A) + ^(r, A)r2(A) dC(A)}. (5.13) 

We have the Pnrseval identity 

oo 

J\gir)\^dr^^\\T\\^. (5.14) 



6. Computing the Measures 

In the sections to follow we will apply the results of section 5 to our particular 
eigenvalue equation, given by (|2.5|l . Thus, for the function Q{r) in H5.1|l we will 
have the expression 

Q{r)^q{r) re 0,oo. 

In this section we will calculate the measures given in H5.9a(l . (|5.9b(l and H5.9c|l . 

Lemma 6.1 Let M^(l), M^{l),'il^'^{r,l), V^(r,0 [l e <C+ = {I : 3ml > 0}) be 
the quantities defined in section 5 for equation (|2.5|l (we use the superscript m to 
emphasize their dependence on m £ ZJ. Let jm(r, /) Wm{r,l) he the solutions of 
(|2.5(l dehned respectively in lemma m\ and by H3.17|l . Then, 

^o"(0 = -^mi^, M::{l)^-^^^y (6.1) 

Jm(RJ) Wm[R,l) 

and 

ro{r,l)^^^ i:Z{r,l)^^^^^^ (6.2) 

Proof. It is easy to note that jm{r, I) and ym{r, I) defined in lemma ITTI are linearly 
independent solutions of (|2.5(l , thus any other solution will be a linear combination 
of these two. Then 9{r, I) and (j>{r, I) defined in section 5 can be represented as 

i0{r, I) = a{l)j„r{r, I) + [3{l)y^n{r, I), 
\ipir,l) = l{l)jm{r,l) + 5{l)ym{r,l), 

for some coefficients a, (3, 7, 5. Let 

A(/) = UR, l)yM - JmiR, i)yM 0- (6.4) 
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Being the Wronskian of two linearly independent solutions A(/) is non-zero. Using 
the boundary conditions (|5.2(l it is easy to find that 

From H6.3|l and (|5.4(l it follows that 

{aj™(i, /) + Py,nit, I)} COST + {af„^{t, I) + (3y'^{t, I)} sinr 



M(0 = - 



{7im(i, + 5ym{t, I)} COST + {-ij'^{t, I) + (5y;„(i, Z)} sinT 



For t ^ the formulas H3.1|) . H3.2a|l and (j3.2b|) tell us that the term ym{t,l) will 
dominate ^-^id jm{t, I), while ym(t, Z) will dominate ym{t, I)- Then, 

Mo™(0 = lhnM(0 = -f. 

By applying (|6.5|l we obtain the first equality in (|6.1() . To get the first equality in 
(|6.2() use definition H5.5|l of tpj{^{r,l) and equalities H6.3|l to H6.5|l . 

The quantities M^{1) and ip'^{r, I) are calculated in exactly the same way. One 
needs to express ip{r, I) and 9{r, I) in terms of Wmir, I) and Xm{r, I), then put t — s- cx) 
in (|5.4|) and use the properties of Wm(r, a;m(r, Z), and their derivatives shown in 
lemma 1^31 The lemma is proved. □ 

The next lemma will study the properties of the functions M™(Z), M™(Z), and 
M™(Z) — M™{1). Here we will set some notation and write some formulas which will 
be used in the lemma. By A we will denote a real variable. Let Wmir, A), km{r, A), 
am{r,\), and bm{r,X) be the functions defined by H3.17|l . H4.5|) and H4.11|l . The 
following equalities hold, 

= Ae-W2j^^(^)^ < argz < n/2 

ITT 

and 

H^l'>{z) = J„i{z) + iY„i{z), -7r<argz<7r, (6.6) 
(formulas (9.6.4) and (9.1.3) from.-) We deduce 

Wmir,X)^—e-'"'/^k^{r,X), X < d\ (6.7a) 

ITT 

and 

Wm,{r., A) = am(r, A) + ib.,n{r, A) A > d^. (6.7b) 



Lemma 6.2 M^{1) is meromorphic across all the complex plane, while M^{1) 
{I £ C+) extends continuously to the real axis. M™(A) — M™(A) is real or inGnite 
for X < (P, and it has a finite number of zeros on the real axis, all in the interval 
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Proof. As stated in lemma im jm{R,l) and j'j^{R,l) will be analytic functions of 
I G C. Then M^{1), being obtained as their ratio, will be a meromorphic function. 

The function Wmif, I) was defined in lemma lX^ It was proved there that Wmir^ I) 
is defined and analytic for all I G {—tt/2 < arg(Z — (P) < 37r/2}. In particular 
Wm{r, X) is defined for all real X ^ cP. To prove that M^{1) extends continuously 
to the real axis, it suffices to show that its denominator, WmiryX), is not zero for 
A e R\{(P} and that lim,^rf2 M^{1) exists. 

If A < , then Wm {r, A) ^ because of H4.5|l and (|6.7a(l , since the Bessel function 
Krnis) takes real strictly positive values for s > (as it follows from section (9.6.1) 
of. II) If A > d^, then w„,.(r. A) ^ because of (jTTTIl and j6Jb|l . since for s > the 
Bessel functions Jm{s) and Yjn{s) take real values and cannot be zero at the same 
time. 

Show that lim;_^£;2 M™{1) exists. If m > 0, then according to (9.1.8) and (9.1.9) 
from, iiwe have that for z G C, z ^ 

H^H^) - - l)!(^/2)"", for m > 1, 

and 

- (-2i/7r)lnz, for to = 0. 
We can extend these to m < 0, by using (|6.t)|l together with 

J-rniz) = (-1)™ J,„(z), r_„(z) - (-l)"y„,(z), TO G Z 

(formula (9.1.5) from.'^ One can obtain the behavior of the derivative of Hm\z) 
as z ^ by formally differentiating the above. Set z = R\/l — d? . Deduce that for 
all m (zl, 

lim M^{1) = H^l/2_ ( 

Note that M™(A) is real or infinite for X < cP, being the quotient of —j'^{R, A) 
and jm{R, A) both of which are real according to lemma IXTl We have that M™(A) 
is real for X < cf, that follows from (|6.7a() and by using again the properties of the 
function X„(s). Then, Mlp{X) - M^'(A) is real or infinite for A < d^. 

Let us show the last part of this lemma, the fact that A/™ (A) — M™(A) finite 
number of zeros on the real axis, all in the interval (0, (P]. Let first prove that this 
function can have no zeros for A < 0. From Ht).7a|l we have 

M^-(\) - M'^(X) = ^-M^ , fcml^.A) ^ A»(i?,A) 

jUR,^) kMX) j,„(i?,A)fc„(i?,A)' 

where 

D^{r, A) = j^„(r, X)k„i{r, A) - j„i(r, X)k'„^{r, A). 

Assume for some A, Mo"(A) - M™(A) = 0. Then D^{R, A) = 0. We will prove 
that is false. First note that Dm(r, A) = Dm{R,X) for r > R. Indeed, jm(r. A) and 
km{r,X) will satisfy (|2.5|l (for jmi^jX) this follows by its definition, for kra{r,X) it 
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follows from the observation that according to (j^.Tafl , km (r, A) is a constant times 
Wm {r, A) which, as defined by 1)3. f 7(1 , is a solution of ((2.5|l for r > R). Then it is 
easy to check that 

D'mir, A) = j::,{r, A)fc„,(r, A) - i,„(r, A)C(^ A) = 0, 

so Dm.{r,X) is constant in r. Second, note that for r large enough Dm{r,X) < 0. 
Indeed, for r > 0, jm(r. A) > 0, A) > by lemma IX^ and fcm(r, A) > by the 

properties of Bessel functions. Also, for r sufficiently large we have k'„^{r,X) < 0, 
this is a consequence of (|4.8|) . We infer that for r large enough Dm{r,X) < 0. 
These two observations imply D,n{R,X) < 0, therefore Dra{R,X) is non-zero, and 
Mo-(A)-M™(A)^0. 

Now let < A < c?^. We showed that Mq"(/) is meromorphic on the whole 
complex plane. The function Wm{r, I), as defined by ((3.17|l . is analytic on {— 7r/2 < 
arg(/ — cf) < 3tt/2}, in particular it is analytic on the set of / such that ^Rel < (P. 
Then M'^{1) = —w'j^{R, I) / Wm{R, I) is meromorphic on the same region, and so 
is M™{1) — M^{1). Therefore, it can have only a discrete number of zeros on the 
interval (0, d^). 

Assuming that the number of zeros is infinite, their only possible accumulation 
point is A = d^. So, there will exist a sequence A„ < (P, n > 1, with limA„ = (P 
and M^{Xn) - A/;:,'(A„) = for all n>l. According to formulas (9.1.3), (9.1.10) 
and (9.1.11) from, -'^ the Hankel hmction H'^^\z) will have the representation 



if«(z) = /i(z) + /2(z)lnz, 

with fi{z) and f2{z) meromorphic functions of z e C. Then, if we recall definition 
(|3.17|) of Wm{r, I), one can calculate that we will have the representation 

Mo-(o-AC(o- ^^i;it^i;!r^ 

with gi,g2,93 and meromorphic functions on the whole complex plane, and 
t = € C. We infer 

gi{tn) + g2{tn) lnt„ = 0, 



for all n > 1, where t„ = \/ A„ — d? . Two cases are possible. If 52(^71) is zero for 
infinitely many n, this implies gi(tn) = at the same points. Then the meromorphic 
functions gi{t) and (72 (i) are identically zero, and so is M^{1) — M^{1), obtaining 
a contradiction. Otherwise, if (72 (in) is zero only for finitely many n, we can write 

lni„ = — -, n>no. 

.92 (in) 

On the right-hand side we have a meromorphic function. As n ^ cx3 we have i„ —>■ 0, 
thus, for some integer p, there follows 

ln|t„| = 0(|i„n. 
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That is clearly impossible. The contradiction shows that Af™(A) — M™(A) can have 
only finitely many zeros on [0, (f). Another possible zero could be at X — (P. 

Lastly, show that Af™(A) — Af™(A) is never zero on {cf, oo). According to 1)6.1(1 
and iHjEI), 



and so, 



Mo™(A) - AC(A) 



{R, X)+ib,n{R,X)' 

jm 0,m ^ ^^m 



(6.9) 



(the arguments (i?. A) were omitted for simplicity) . If we assume that for some 
X>(P, M^'{X) - M™(A) = 0, this implies 

j^(i?, A)a„(i?, A) -a^{R, X)j„r{R,X) = 0, 
j'^iR, X)bM A) - 6'„,(i?, X)j,n{R, A) = 

(since these quantities are real, as it follows from lemma rOl and H4.11|l l. The 
numbers jm{R,X) and j!^{R,X) cannot be both zero, since jm{'r,X) is a non-zero 
solution of the second order differential equation 12.5|l . In this case the vectors 
{ajn{R,X),a'j^{R,X)) and (&m(i?. A), 6^(i?, A)) must be linearly dependent. That 
cannot be since the functions am(r. A) and 6m(r, A) (r > R) are also solutions of 
(ESI, and they are linearly independent. Thus M^(A) - M™(A) ^ 0. This finishes 
the lemma. □ 



Theorem 6.1 Let ^,„, rjm and Cm be the functions defined by (|5.9a|) . (|5.9b|) and 
(|5.9c|l . There exists a non-decreasing function Xm : M ^ R such that the following 
measures are equal 

d^m{X) = jni{R, Xf dXm{X), 

dVmW = -jM A)jm(i?, A) rfXm(A), (6.10) 
dCmiX)^j'^iR,XfdXmiX). 

The function Xm is identically zero for X G (— oo, 0], is piecewise constant for X G 
(0, d^) where it has a finite number of discontinuities, and is continuous for X G 

(d2,0o). 

Proof. Denote 
According to (|5.9a() . for any Aq < Ai real numbers 

Ai 

Cm(Ai) - ^m(Ao) lim 3mM{s -i-iS)ds. 

5^0+ J 
An 
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In particular, if M{1) extends continuously to the interval [Ao,Ai], then by using 
Lebesgue's theorem of dominant convergence it is easy to show that 

Ai 

em(Ai)-U(Ao) - JjmM{s)ds. (6.11) 

Ao 

The same kind of reasoning clearly holds for rjm and Cm . 

As it follows from lemma 1^!^ if A < then M™(A) — M™(A) is real or infinite, 
and non-zero. By applying H6.11|l we find that for any Ai < A2 < 0, ^m(Ai) — 
^„(Ao) = 0. From (|5.9a|) we have that ^m(O) = 0, and thus Cm (A) = for all A < 0. 
In the same fashion one obtains that for A < 0, ?7m(A) = and Cm (A) — 0. Set 
Xm(A) = for A < 0, and (jOTjl will hold. 

Let A'" < A™ < • • • < A™ be the points in the interval (0, (P] where, ac- 
cording to lemma lO M|5"(A) - M™(A) = 0. We can use lemma lO and (|SjT|) 
to deduce that Cmj Vm and Cm are constant on each of the intervals making up 
(0, (i^]\{A™, Af*, . . . , Ap^}. Set Xm to be constant on each of these intervals. At 
each of the points A™, A™, . . . , A™ the functions Cm, Vm and Cm could have a jump. 
To show (I6.10() on (0, (P] we need to find a relationship between the jumps of these 
functions. We will return to this shortly. 

The remaining case, A > c?^, is treated similarly. Equation 1)6. 9|l in the proof 
of lemma 1^21 gives an expression for M™(A) — M™(A) on this interval (with the 
notation from (|4.1Hl 'l. By using the fact that the quantities a„i, and jm together 
with their derivatives are real, we can calculate 

3mM(A) = Jm|- ^ Jm(am + »b ) ^ 1 

Use ()4.12() to simplify the numerator of this fraction. Apply (|6.11() . We get that for 
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any d"^ < Aq < Ai, 



U{\i)-U{Xo) - 2ic„.(A)2+d„(A)2^^' 

Ao 

where Cm(A) and dm (A) are defined by H4.13a|) and l|4.13bp . So we have 



TT j„,{R,\ydX 



2c„(A)2+d„(A)2- 
We want (|6.10|l to hold. Define Xmi\) for \ > such that 



2c™(A)2+dm(A)2' 

then the first of the three identities (|6.10() is valid. It is easy to repeat the same 
calculation for rjm and Cm and show that for A > d^ the other two identities in 
hold. 

Now we will return to what is the longest part of the proof, the study of what 
happens at the points A G (0,c?^] where Af™(A) — Af™(A) = 0. Let Aq be such 
a point. An immediate observation is that AfQ"(Ao) is finite and jm{R,^o) 7^ 0. 
Indeed, it was proved in lemma [^7^ that M™(A) is finite for A real. So, Af™(Ao) 
which equals M™(Ao) is finite. Then, since M™(Ao) = — Ao)/jm(-R, Ao) and 
because j'„(-R, Aq) cannot become zero simultaneously with jm{R, Aq) {jmi^, I) is a 
solution of I2.5|l '). we deduce jm{R, Aq) 7^ 0. 

With this observation in hand, in order to prove that l|6.10|) holds at Aq one 
needs to show that 

dVmiXo) = Mo"(Ao) d^iXo), dCm(Ao) = Mo"(Ao)' dU{Xo), (6.12) 

and then define dXm(Ao) = di]^{Xo)/jm{R, Aq)^- 

Let vq be the jump of at Aq. Recall, was defined by (|5.9a|l . so, 

Ao+6 

f 1 

rn = lim lim / —Urn -, ^ z rr ds. 

" e^o+s^a+ J M^'{s + i5)- M^{s + i5) 

Ao-e 

By using Ij5.9bp . the first equality in lt).12|l can be written as 

Ao+e 

/A//"™ i s -I- lA") 

Ao — £ 

M^{1) is analytic around Aq. Therefore, in a neighborhood of Aq 

A/^"(s + i6) = Mo"(Ao) + {A/o"(s) " ^^™(Ao)} + ^SH{s + iS), 
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with H a continuous function around Aq. Substitute this above. By using the fact 
that M™(Ao) is real and the definition of tq, we get the equivalent equality 



Mo™(Ao)ro 



lini lini 



Ao-£ 



{Mo"'(g) - M^jXo)} + zSHjs + tS) 

= Mo™(Ao)ro, (6.13) 



so we have to prove that the limit of the integral on the left-hand side of H6.13|l is 
zero. To do this, we will need additional information about — M™. 

Let I ~ s + iS {6 > 0). Recall that formulas (|5.6() and H5.8|l hold, where ip™ and 
■0™ are defined in (|5.5|l and (|5.7() and an expression for them is given by (|6.2() . Add 
equalities 15.6|l and 15.8|l . Obtain 



.3m{Miril)-MZm 



jrn{r,l) 



dr - 



WmirJ) 



W,niR,l) 



dr. 



This equality gives us two things. First, that 3m{M™(/) - M"(Z)} > 0, therefore, 

1 



Mo"(0-m™(0 



> 0. 



(6.14) 



Second, for I close to Aq, the quantity (5^^3m{Af™(Z) — M™(/)} is bounded from 
below by a strictly positive number, say w^^. Then (5^^ |M™(Z) — M™(Z)| is bounded 
below by the same number and therefore. 



|Mo™(0-Af™( 
The integral in ()6.13fl can be written as 



(6.15) 



Ao+e 



'3m 



Ao-e 



M^js) - M^"(Ao) 
M^{s + id)- M^{s + iS) 



ds 



Ao+e 



-am 



iSHis + iS) 



Ao — e 



ds. (6.16) 



Denote q = sup|3_;)^^|<j |Mg"(s) — A/™(Ao)|. Since Af™ is continuous in a neigh- 
borhood of Ao, we will have as e — > 0. Let us estimate the first integral from 
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— Jm 



Ao-e 



Ao+e 



Ao+e 



ds < 



Ao— e 



|M„"(s)-A.C(Ao)| 



Jm 



3m 



An — e 



Ao+e 



3m 



Ao— e 
Ao+e 



Ao— e 









f iS) 


- M™(s + z(5) 






1 






-MJS{s + iS) 






1 




f i<5) 


- A/™(s + i6) 






1 










(Ao + 





ds 



ds 



ds 



In deriving this we used the fact that A/™(s) — AfQ"(Ao) is real, together with (|6.14() 
and H5.9a|l . Clearly as £ ^ 0, the integral goes to zero. 

Now estimate the second integral in the sum (|6.16|) . Let i/oo be an upper bound 
of \H{s + id)\ for I — s + iS in a neighborhood of Aq. Inequality ()6.15f) implies that 
as e — > 0, 



Ao+e 



-3m 



iSH{s + id) 



Ao-e 



M^{s + i5)-M^{s + iS) 



ds 



Ao+e 



< 



Hca^oo ds = 2eHaoUJa 



This proves the first equality in H6.12|l . 

Let us prove the second equality in H6.12(l . Recall that is given by H5.9c|l . The 
above approach does not apply immediately, since unlike M™{1) (the numerator in 
(l5.9b|l V the function M^{l)M^{l) (the numerator in 15.9c|l 'l will not be analytic 
around Aq if Aq = (since as seen from lemma lH.41 ^^(r, is not defined in a 
neighborhood oi I — d^). The idea is then to use the equality 

xy 



X — y X — y 
to write the jump of Cm at Aq as 

Ao+e Ao+e 

/M"^(s + iS) f 
-3m—— ^^r:^ 7Trds-\- lim lim / 3m (s iS) ds . 

Ao-e Ao-e 

The limit of the second integral in the sum will be zero, since as argued above, 
M™ will be finite at Aq (and therefore, around Aq) and thus the quantity inside the 
integral is bounded. For the first integral in the sum we can proceed in the same 
way we calculated the jump of rjm- This finishes the proof of (|6.10|l . 



26 Wave Propagation in a 3-D Optical Waveguide 



Finally, we need to justify the claim that Xm is a non-decreasing function. From 
we have 

dUW+dCm{\) - {j™(i?,A)2 +j;^(i?,A)2}dx™(A). 

The left-hand side of this is non-negative measure, since by theory and Cm are 
non-decreasing. The number jm{R, A)^ + j'^iR, A)^ is strictly positive, as jm{R, A) 
and Jm(i?, A) cannot be both zero, jm{r,X) being a non-zero solution to the sec- 
ond order differential equation (|2.5|l . Then we get that dxmW is a non- negative 
measure, which shows that Xm is non-decreasing. □ 

Corollary 6.1 Let A e (0, rf^] be a discontinuity point for Xm- Then, if A < d^, the 
following hold: 



3m 

(i?,A) fc™(i?,A)' 



(6.17a) 



and 



While for \ = d^ 



j,n(r,A)^ f"i^'^| fc,„(r,A), r > R. (6.17b) 
k,n(R, A) 



j;„(fl,A) _ |m|-l/2 
j™(i?,A) i? ' 



(6.18a) 



and 



Jn particular, for A G (0, d^] a discontinuity point of Xm, jmir, A) decays exponen- 
tially as r —> oo if A < d^, and jmij", A) ~ as r ^ oo if A = d^. 

Proof. As it follows from the proof of Theorem 16.11 for such a A we will have 
Mo"(Ao) - M"(Ao) = 0. With the help JHIIl), I^Tajl and (ESJ, we deduce (|CT7al) 
and (|6.18a|l . Notice that these equalities are exactly the conditions H4.6|l and (|4.9I) . 
Then H6.17bp and (|6.18b|) follow from (gTl and □ 

7. Computing the Transform 

Denote by L^(xm) the space of all functions G : M ^ C such that 

|G(A)|2dx™(A) <oo, 

where x„i is the non-decreasing function defined in Theorem 16. II 
Theorem 7.1 Let g e L'^{0,oo). The integral 

oo 

G™(A) = ( 3,,,{r,\)g{r)dr (7.1) 
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is convergent in LF'{xm), in the sense that there exists Gm € L'^iXm) such that 
— > Gm in L'^{xm) as c — > and oo, where 

d 

GmW = Jjrn{r,X)g{r)dr, < c < d < ^. (7.2) 

C 

The equality 

oo 

<?M-^ I j^{r,X)Grn{X)dxm{X) (7.3) 

— OO 

holds, in the sense that g^^ g in L^{0, oc) as r ^ — oo, a oo, where 

a 

g^"{r) J jrn{r, A)G™(A) dx,n{X), ~oo < a <t <oo. (7.4) 

r 

We have the Parseval identity 

oo oo 

J |3(r)|2 dr = i 1 |G™(A)|2 dx^W- (7.5) 

-oo 

Proof. We will apply TheoremlO First note that if T = (ri,r2) with Ti, r2 : 
R C, then because of Ijti.lOII . the norm || • || as defined in (|5.11|) can be written 
in the form 

oo 

l|r|p= [ \jm{RA)riiX)-jURA)r2iX)\'dxmix). (7.6) 



Recall the identities (|6.3|l to (|6.5|) for expressing d{r, A) and (p{r, A) in terms of 
jm{r, A) and ymir, A). Note that from (|6.4|) and (|6.5|) we get the following equalities 
involving the coefficients a,f3,"f, and d 

a{X)jM X) ~ i{X)3'„XR, X) = 1, (7.7) 

and 

P{X)jm{RA) - 5{X)j'm{R,X) = 0. (7.8) 
Then, for < c < d < oo, the hmctions V\'^ and defined by (|5.12|) become 



rf (A) = j {a{X)3m{r,X)+mVnr{r. X)]g{r)dr, 

C 

d 

n\X) = I {7(A)j™(r, A) + 5{X)yra{r, X)}g{r) dr. 
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We denoted T^'' = (rf ,rf ). Let 

d 

f f (A) = / a(A)j™(r, X)g{r) dr, f f (A) = / 7(A)j™(r, A)g(r) dr. (7.9) 



and set T'^'^ = {rf,Tf). We have 

||pC<i_fcd|| ^ Q_ (7.10) 

That follows from (|7.6|l and l|7.8(l . Theorem 15.11 guarantees the existence of F = 
(FijTa) e L^{p) such that 

I |r - r^'^l I ^ as c ^ 0, ^ oo. 

Then, equality H7.10|l says that we have ||r — f'^''|| ^ as c ^ 0, d ^ oo. But, 
according to l|7.6|l . 

oo 

l|r-^1P= J \{Ur,x)t,{x)^jUR,x)t,{x)}- 

— oo 

{j,n{R, A)f f (A) - jUR, A)f ^^(A)}|2 dxm. (7.11) 

Note that 

d 

A)f f (A) - jUR, ^)rfiX) = / j;„(r, A)g(r) dr. 



That follows from (fT7|l and ((7^ . Thus, if we denote 

G„(A) =j™(i?,A)ri(A)-i:„(i?,A)r2(A), AeR, (7.12) 
we obtain from H7.11() that 

\Gm{X) - G^nWl^ dXm ^ as C ^ 0, d ^ OO 

(G^ was defined by H7.2(l 1. This shows the first part of Theorem 17. II 

Next we need to show that representation H7.3|l holds. It suffices to prove that 
g'^"' as defined by H5.13|l in Theorem 15. II is the same as g'^'^ defined in (|7.4() . And 
they are. To check this one needs to start with 5'^'^ as given in H5.13|l . substitute 
0{r, A) and (p{r, A) from H6.3|l . use H6.10|l to express ^m, rjm and Cm in terms of Xm, 
use the equalities 1)7. 7|) and H7.8|) . and finally use definition (|7.12() for Gm(A). 

Lastly, the Parseval identity l|73|l follows from (jHHH), ((7^ and (fTT^ . The 
theorem is proved. □ 
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Theorem 7.2 Let g e L^(0, oo). Let Xm be the non-decreasing function defined 
in tlieorem^i^ Let < A™ < • • • < AJf^^ < {Pm > 0) be the points where Xm 
is discontinuous. Let r™,...,rp^ be the corresponding jumps. Let amir, X) and 
bm{r,X) be the functions defined by (|4.1HI . and Cm{X) and dm{X) be defined by 
(lOl^ and (l4l3hll . Then, 



-1 

\m\2 







rT = ^i 3rn{r,X^Ydr\ ,k=l,...,Pm, (7.13) 



'^^'"W = o -TTTy— TTTW' Ae(d2,cx)). (7.14) 



vr dX 
2c™(A)2 + rf„(A) 

We Jiave tie representation 



p oo 



k=i 



Proof. In theorem 16.11 we proved the existence of the function Xm and along the 
way we found its continuous part, that is, equahty (|7.14|l . We need to find its 
discrete part, that is, the value of all the jumps of Xm- Then H7.15|l will follow by 
applying (|7.3() . 

Let us notice the following observation. If Ai ^ Xi, and w\(r) and wi(r) satisfy 
(12. 5() with A = Ai, and A = A2 respectively, then for any < c < 0? < 00 



d 

iL..' , 



wi {r)w2 (r) dr = -{Xi - X2) [w[ {r)w2 (r) - wi (r) (r)] 



(7.16) 



Indeed, we can write 

d d 



Tn?" — 1 /4 ' 

'w'l{r)w2{r) dr ^ - I IXi- q{r) \wi(r)w2{r) dr. 



and 



d d 

f ( - 1/4", 
Wi{r)w2{r) dr ^ - Wi{r)< X2 - q{r) \ w2{r) dr. 



If we integrate by parts the left-hand sides of these two equalities, and then subtract 
from first the second, we get exactly (|7.16|l . 

Before we prove (|7.13|) . recall the behavior of jm(r, A) as r — > oo. For A < a 
discontinuity point of Xm it is described in corollary 16. II while for X > d^ see the 
discussion in section 4. 

Let Ao < d^ be one of the discontinuity points of Xm- Let rg be the corresponding 
jump. We will consider two cases: when im('', Aq) is square integrable, and when 
it is not. The first case splits into two subcases: we can have either Aq < d"^, or 
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Aq = (P with |r7i| > 2. The second case happens for Aq = (P and \m\ G {0, 1}. Let 
us start with the first case. 

Denote g{r) — jm(r, Aq). Since g{r) is square integrable, we can apply theorem 
17.11 for this particular function. We will show that the corresponding Gm(A) as 
defined by H7.1|l is such that 



for all A such that dx,n(A) ^ 0. Then (|7.13(l will follow promptly; one needs to apply 
the Parseval identity (|7.5(l and notice that since tq is the jump of x,„ at Ao, then 
(iXm(Ao) = roS{X — Aq) with S being Dirac's function. 

Consider first the subcase Aq < cP. That H7.17|) is true for A = Aq follows 
immediately from (|7.1(l . Assume now A 7^ Aq. Let us compute G"l{\) for < c < 
d < CX3 as defined in (|7.2|l . Use (|7.16|l with wi[r) = jm{r, Aq), W2{r) = jm{r, A). Put 
c ^ and d ^ 00. From (|3.1(l we deduce that 



Also, both jl^{d, Ao)jm(d, A) and jmid, Ao)j^(c?, A) go to zero as d — > 00 if dxmi^) 7^ 
0, since on one hand, jmif, Aq) and its derivative decrease exponentially, and on the 
other hand, j'^^ir, A) and its derivative either decrease exponentially for A < d^, or 
behave like a power of r for A = d^, or are bounded for A > d^. In any case we get 
that G^(A) as c and d ^ c», so G„(A) = 0. 

The subcase Ao — d^ with |m| > 2 follows in the same way. The statement 
that for A Ao and dx„(A) ^ both j^(d, Ao) jm(d, A) and jm{.d,\o)j'^{d,X) go 
to zero as d — + cxD is argued in a little different way. We have that jm{r, Ao) and 
its derivative decay as a negative power of r for r — > 00, while jmir, X) and its 
derivative either decay exponentially for A < Ao, or stay bounded for A > Aq. But 
the conclusion is the same, G^f (A) — > as c — * and d ^ 00, and thus (|7.17() holds 
in this case too. 

Now consider the second case, Aq — d^ and |m| G {0, 1}, when, as we remarked 
above, jm{i", Ao) is not square integrable. For < c < d < 00 define 



This function will be square integrable. Let Gm(A) be the corresponding transform 
of g{r) as defined by H7.1|) . Apply the Parseval identity H7.5|) . Get 




(7.17) 



im(c, Ao)im(c, A) - j™(c, Ao)j™(c, A) ^ as c ^ 0. 




jm(7', Ao), if c < r < d, 
0, otherwise. 



d 



00 




c 



— 00 




{Ao} 
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From (j7.1|l we obtain that 

rd 



Then we can write 



G„,(Ao)= / \jmir,\o)\'dr. 



(r, Ao)pdr^ > ro. 

By putting c — > 0, d ^ oo and noticing that rg > 0, being a jump of the non- 
decreasing function x,„, we deduce tq = 0. □ 

With this theorem, coroUarv lfi. II and the classification obtained in section 4, we 
can characterize completely the points of discontinuity of Xm- 

Corollary 7.1 Let A S K. Then A is a discontinuity point of Xm if and only if 
jm(-, A) G i^(0, oo), and if and only ifO<X<(f and (|6.17a|l holds, or X — d"^ and 
(|(j.l8a|l iio7ds. 



8. Finding Green's Function 

In this section we will show that under certain conditions, the solution of the 
Helmholtz equation which in the cylindrical coordinate system (r, z?, z) is 

written as H2.1|l . is unique. We will find a representation for it in terms of the 
source /(r, z), the eigenfunctions jm{r, X) of equation l|2.5|l satisfying lemma IXTl 
and the measure dXm(A) defined in theorem l6.II Before proving this result we will 
need one lemma. 



(8.1) 



Lemma 8.1 Let u E C^(K'^). Then u can be written as 



For each z G M, the function r Um{r, z) is in C^[0, oo) and 
d{^/rum{r, z)} djm{r,X) 



lim 



jrn{r, X)- 



dr 



dr 



-{^/rUm{r,z)} 



= 0. 



(8.2) 



Proof. Equality (|8.1|) is nothing but the Fourier series of the function •& u{r, -d, z). 
The smoothness of the obtained Um follows from the formula for the Fourier coeffi- 
cients. 



2tt 



1 



Umir,z) = I uir,t,z)e-'"''dt. 
Zn 



Let us prove 

(IH3- Write it as 



lim 



Vrjni{r, X) — ^^-^ h Wm(r, z) 



im{r,X) djrn{r,X) 



dr 



0. 
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The first term in the sum clearly goes to zero as r — s- 0. For the second term, by 
applying lemma mi to the expression in parentheses we get 



lim / jm(r,A) djm{r,X) 



dr 



1 



lim f ^rl™l -<\m\ + ^V'™') = - lim |m|rl™l = 0, 



for all TO e Z. □ 

These will be the conditions on the solution u for H2.ll) which will guarantee its 
uniqueness. First, we will assume that the source / is continuous and with compact 
support. Second, we will impose the condition that u G C^(K-^). Third, suppose 
that for all TO e Z, z g M the following equality holds 



lim 

r — >oo 



Jm{r, A) {Vru„i{r, z)} 



dr 



dr 



= 0, 



(8.3) 



with the functions Um{r, z) defined by H8.ll) . 

Denote by Um{X, z) the transform of the function r y/rum{r, z) given by (|7.1() . 



Um{\ -z) = / ira{p, X)\/rUm{r, z) dp, 



(8.4) 



The fourth requirement is the radiation condition 



lim 

1^1 — >oo 



0, for A < Pnl with dx-m ( A) 7^ 0, 



lim Um{\ z) = 0, 



for A > fc^ri-Q. 



(8.5) 

These conditions are physically motivated. First condition says that the source 
must be finite. Equation H8.3|) signifies a fast decay of the electromagnetic field 
intensity as r — > c». And the radiation condition H8.5|l means that the energy 
going to 2; = cxD can be separated in two parts. First part is oscillatory, and it 
goes to infinity, and cannot not come from infinity, while the second part is rapidly 
decaying. 



Theorem 8.1 With the above assumptions, the solution of H2.1|l can be repre- 
sented as 

00 00 27r 

u{r,^,z)= J J jG{r,p;^,t;z,C)fip,tX)pdtdpdC, (8.6) 

-00 
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where 



1 1 ^r°pi\^-(\\/k^nl-X 

^ / ^ e-("-*)j^(p,A)j^(r,A)dXr„(A), 



< r,p; < < 27r; z,C e K, (8.7) 
and Xm is the non- decreasing function defined in Theorem \6.1\ 
Proof. The function f{p, t?, z) can be decomposed as 

^e'™Vm(p,^), (8.8) 

with 

1 



fm{p,z)^^J f{p,t,z)e-'^'dt. (8.9) 



Look for u{p,'d,z) in the form (|8.1|l . By plugging (|8.1() and ()8.8|) into H2.1|) we 
deduce that for each m, Um{p, z) needs to satisfy the equation 

for aU TO e Z. Let A G M be such that dxm{^) 7^ 0. Let Um{^, z) be the transform 
of u„i{p, z) given by H8.4|l . and let 

C30 

F™(A,z)= jm{P,>^)y/pfm{p,z)dp (8.11) 



be the transform of fm. Multiply (|8.10|l on both sides by y/pjm{p, A) and integrate 
from to CX3. Obtain 



d^Urr 



oo 

, /". / .^ ^ d f dum{p,z) \ 



dp+ 



jniiP, A) (^'^n{pf - ^u,n{p, z) dp = i^„i. 



Use integration by parts twice for the first integral in the above equation. By 
applying lemma IHTI and equality (|8.3|) we get 

"dz^ J \ dp V ljm(/5,A)j-^u„i(p, A)dp = F„j. 
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Recall that jm{p, A) satisfies H2.5|l with q{r) given by (|2.3|l . Then, 

d^U„ 



+ {k'n'a-X)Um^F^. (8.12) 



9z2 

The solution to (|8.12|) which satisfies H8.5|) is easily found, 

— oo 

or if we use H8.11|l . 

+00 00 



UmiX,z)= I I ^^^—j===-jrnip,X)^/pfni{p,OdpdC. 



-00 



Um{X, z) was defined by (|8.4|l . Using the inversion formula ()7.3|) given in theorem 
17.11 we can recover u,n{r, z). 



00 

/rw„(r, z) = - / jm(?', A)t/m(A, z)dxm(A), 



00 +00 00 



u,n{r,z) = ^— I I I " ^^^ ^j„,ip, \)jm{r, X)pfm{p,0 dpdC dxrn{X)- 
T^^frp J J J 2i^Jk^nl - X 

—00—00 

Now, to get H8.6|l with G(r, p; ??, 2;, Q) given by H8.7(l we need to substitute fm{p, z) 
from H8.9(l . find u{r,-d, z) from H8.1|l and interchange the order of integration so that 
the inner-most integral is the one in respect to A. □ 

The theorem we just proved shows that the electromagnetic field generated by 
the source / can be decomposed in two parts: the guided part, which is a sum of 
guided modes decaying in r either exponentially or as a power of r, and a radiation 
part, which is obtained by summing in m and integrating in A. For each m G Z the 
set of guided modes is finite, as it was shown in lemma 1^?^ The next theorem will 
prove a stronger result. 

Theorem 8.2 The total number of guided modes (in all m G Z) is finite. 

Proof. We just need to show that for |m| large enough there are no more guided 
modes. A mode jm(r. A) is guided, if 

AC(A)-M™(A)=0. 

All A for which this equality happens are in (0, cZ^], as proved in lemma By 
using (|6.1|l . (|6.7a|l and H6.8|l we can write the this equality as 

j'^{R,X) ^ k'„^{R,X) ^ _^ ^2 
(i?,A) fc™(i?,A)' 
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or 



jm{R, A) 



\m\ - 1/2 
R 



, if A = c^^ 



where k{r, A) is given by 14.5|l . 

The left-hand side of these equalities is strictly positive for \m\ large, as it follows 
from lemma IX^ Thus, the second of these equalities is not possible. To show that 
the first one cannot happen, it suffices to prove that k'„^{R, A) < for |m| > 1, since 
we know that km{r, A) > 0, for all r > 0. 

The function r — > km,{r, A) will satisfy the equation 



which implies that k'^^{r, A) > for r > 0, and so, k'^{r, A) is an increasing function 
of r. From this and 14.8|l it follows that A;^(r, A) < for all r > 0. This finishes the 
proof of the theorem. □ 

9. Conclusion 

In this paper we have constructed a framework for analyzing waveguide problems 
which is based on a transform theory. The construction of the transform was more 
difficult, but the final form relatively similar, to the 2-D case.'^The primary tool in 
obtaining the transform was the theory of self-adjoint singular eigenvalue problems. 

This paper completes the study of the wave propagation in a infinite cylindrical 
waveguide. We obtained a Green's function valid for every choice of the index of 
refraction of the core with cylindrical symmetry. In particular, it is enough to solve 
(|6.17a|l for A S (0, d?] and the differential equation (|2.5|l in the core region in order 
to obtain the corresponding Green's function. 

The obtained formula for Green's functions is very amenable to computation. 
In a future article we will calculate explicitly Green's function in the cases of a 
step-index fiber and a coaxial waveguide and will display numerical results. 
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